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Abstract 

We conducted astrometric VLBI observations of water-vapor maser emission in the massive star forming 
region IRAS 21379+5106 to measure the annual parallax and proper motion, using VERA. The annual 
parallax was measured to be 0.262 + 0.031 mas corresponding to a trigonometric distance of 3.82lg' 5 ^ kpc. 
The proper motion was (/ u a cos6,fis ) = (-2.74 + 0.08,-2.87 + 0.18) mas yr _1 . Using this result, the Galactic 
rotational velocity was estimated to be Vg = 218 ± 19 km s _1 at the Galactocentric distance R = 9.22 ± 0.43 
kpc, when we adopted the Galactic constants Rq = 8.05 + 0.45 kpc and Vo = 238 + 14 km s _1 . With 
newly determined distance, the bolometric luminosity of the central young stellar object was re-evaluated 
to (2.15 + 0.54) x 10 3 Lq, which corresponds to spectral type of B2-B3. Maser features were found to be 
distributed along a straight line from south-west to north-east. In addition, a vector map of the internal 
motions constructed from the residual proper motions implies that maser features trace a bipolar flow and 
that it cannot be explained by simple ballistic motion. 

Key words: astrometry, H n regions, Galaxy: kinematics and dynamics, ISM: individual (IRAS 
21379+5106), stars: formation 


1. Introduction 

VERA (VLBI Exploration of Radio Astrometry) is 
a Japanese VLBI system dedicated to exploring the 
three-dimensional Galactic structure, being operated by 
the National Astronomical Observatory of Japan and 
Kagoshima university (Honma et al. 2000). It consists of 4 
stations, each with a 20-m antenna, located at Mizusawa, 
Iriki, Ogasawara, and Ishigaki. In order to achieve high- 
precision astrometry, a dual-beam observation system is 
employed to observe Galactic maser sources and a distant 
QSO as a reference source, simultaneously (Kawaguchi et 
al. 2000). Consequently, real-time phase calibration in¬ 
creases the signal-to-noise ratio and improves positional 
accuracy. 

The “Outer Rotation Curve” (ORC) project is one of 
the main projects carried out with VERA, which aims to 
derive the rotation velocity of the outer part of the Galaxy 
(R > i?oi R denotes Galactocentric distance and subscript 
0 denotes the Sun). Since the inner rotation curve can 
be measured using velocities at the tangent points, it has 
been well determined (Burton & Gordon 1978; Fich et al. 
1989; McClure-Griffiths & Dickey 2007). However, in the 
case of the outer rotation curve, there are no such useful 
points for isolating the terminal velocity. Therefore, the 
distance of each outer Galactic source has to be measured 
directly. As such, in comparison with the inner Galaxy, 


there is much larger uncertainty seen when compiling the 
available outer rotation curves (Dehnen & Binney 1998; 
Sofue et al. 2009). More efforts such as VLBI astrometry, 
are necessary to obtain a more accurate outer rotation 
curve for the Milky Way Galaxy. 

IRAS 21379+5106 is one of the ORC project tar¬ 
gets, which is known to be an ultra-compact H n region 
(Bronfman et al. 1996; Wang et al. 2009). Since ultra 
compact H n regions are thought to be formed around 
massive stars, with a main sequence luminosity of B3 or 
hotter (Churchwell 2002), IRAS 21379+5106 is consid¬ 
ered to be a high-mass star-forming region (Wu et al. 
2006; Shepherd & Churchwell 1996). A far-infrared, three 
color (25, 60 and 100 /iro) image of the area around the 
source is shown in figure 1 , where IRAS 21379+5106 is 
a point source indicated with a green circle. The position 
is (a, 6) = (21 h 39 m 40?8, +51°20 / 35") (J2000.0) in equato¬ 
rial coordinates (Sunada et al. 2007), which corresponds 
to ( l,b ) = (95?29670, —0?93667) in Galactic coordinates. 
The radial velocity with respect to the Local Standard of 
Rest (LSR) is V LS r = -42.27 + 0.22 km s" 1 from 13 CO 
J = (2 — 1) observations (Wang et al. 2009), which we 
use as the LSR velocity of the source throughout this pa¬ 
per. IRAS 21379+5106 is thought to be associated with 
the Perseus arm, considering its position in Galactic co¬ 
ordinates and the LSR velocity (Dame et al. 2001), and 
the heliocentric distance has been estimated to be 5.9 kpc 
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based on kinematice distance (Wang et al. 2009). 

This source has been detected in some molecular lines 
such as CS( J =2-1) (Bronfman et al. 1996), 12 CO(J =1- 
0) (Wouterloot & Brand 1989; Shepherd & Churchwcll 
1996), 13 CO(J =2-1) (Wang et al. 2009), H 2 0 maser 
(Wouterloot et al. 1993; Brand et al. 1994), and NH 3 
(Sunada et al. 2007). Wouterloot & Brand (1989) noted 
that the 12 CO(J =1-0) spectrum showed a red shoul¬ 
der profile, and the 13 CO(J =3-2) spectrum showed two 
peaks, with brighter emission in the blue-shifted compo¬ 
nent (Wang et al. 2009). These line profiles can be in- 
terepreted as bipolar molecular outflows asscociated with 
IRAS 21379+5106. 

As Part II of the ORC project series, we report measure¬ 
ment of the annual parallax and proper motion of IRAS 
21379+5106, observed with VERA, as well as the distri¬ 
bution and kinematics of the maser features. Section 2 de¬ 
scribes details of the astrometric observations with VERA. 
Details of data reduction are summarized in Section 3. 
Section 4 presents our results, including: the annual par¬ 
allax, trigonometric distance, proper motion, and maser 
distribution map. Section 5 discusses the 3-dimensional 
(3-D) motion of IRAS 21379+5106 in the Galactic disk, 
modification to the spectral type of the central young stel¬ 
lar object (YSO), and a bipolar outflow seen in the maser 
distribution. Section 6 reports the summary. 

2. Observation 

VLBI monitoring observations of IRAS 21379+5106 
were carried out in 9 epochs from 2009 November 6 to 
2011 November 5 in K-band with typical observation time 
of 8 hours. The observation frequency was set to 22.232- 
22.248 GHz to detect 6i6 —» 5 23 transition of the H 2 0 
maser, which has a rest frequency of 22.23508 GHz. 

The observations were conducted in dual-beam mode 
(Kawaguchi et al. 2000); one of the two beams, re¬ 
ferred to as the A-beam, was pointed to the main 
target, IRAS 21379+5106, and the other beam, the 
B-beam, was pointed to the phase reference source 
J2137+5101. The tracking center was set to (a, 6) = 
(21 h 39 m 40?55, +51°20'34."0) (J2000.0). The separation 
angle of the sources was 0?52 and both targets were ob¬ 
served simultaneously for phase calibration. The system 
noise temperature was typically 100-200 K throughout 
the observations. In addition to these two sources, we 
observed the bright quasar 3C454.3 with both the beams 
every 80 minutes for fringe finding and band-pass calibra¬ 
tion. 

Other maser sources, IRAS 23385+6053 and IRAS 
22555+6213, and continuum sources, J2339+6010 and 
J2302+6405, were also observed during these observation 
runs for other studies. During each observation run, each 
source was observed for an integration time of 10 min¬ 
utes every 20 minutes. Hereafter we focus only on IRAS 
21379+5106 since the other targets are outside the scope 
of this paper. 

The received signal was output into 16 Intermediate fre¬ 
quency (IF) channels after being filtered with the digital 
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filter unit (Iguchi et al. 2005). The bandwidth of each IF 
was 16 MHz and signals were recorded with 2 bit sampling 
at a rate of 32 MHz. The total recording rate was 1024 
Mbps. The correlation was performed with the Mitaka 
FX correlator, which outputted the auto-correlated and 
cross-correlated spectra with an interval separation of 1 
second. 

One IF was assigned to the A-beam, observing the 
maser line emission. The time-domain data was Fourier- 
transformed into spectra with 512 channels and a channel 
separation of 31.25 kHz, which corresponds to a velocity 
spacing of 0.4213 km s _1 . The other 15 IFs, each with 
a bandwidth of 16 MHz split over 64 channels, were as¬ 
signed to the B-beam, observing continuum emission from 
the reference sources. During data reduction 14 of the 
continuum IFs were used and the 15th IF was abandoned 
due to a difference in correlator settings. 

Observational settings and scheduling were compiled in 
the VEX (VLBI Experiment) format. Basic information 
for all the observations such as epoch, VEX file, Day-of- 
year (DOY) from 2009 January 1, synthesized beam, and 
noise are summarized in Table 1. 

3. Data Reduction 

Data reduction was done with the software AIPS 
(Astronomical Image Processing System), developed by 
NRAO (National Radio Astronomy Observatory), which 
is commonly used in the data reduction of radio astron¬ 
omy. 

The basic procedure of the primary reduction was 
as follows; (1) loading fits files (AIPS task: FITLD), 
(2) sampling bias correction using auto-correlation data 
(ACCOR), (3) a prior amplitude calibration (APCAL), 
(4) correction of delay-tracking, which is necessary to 
replace the tracking model used in Mitaka FX correla¬ 
tor with more accurate tracking model calculated with 
CALC3/MSOLV (Jike et al. 2005; Manabe et al. 1991), 
using a calibration tabic provided by the observatory 
(TBIN), (5) fringe finding of B-beam data (FRING), (6) 
applying the result of (5) to A-beam data (TACOP), 
(7) self calibration using B-beam data (IMAGR, CALIB, 
CLCAL), (8) applying the result of (7) to A-beam data 
(TACOP), (9) applying the calibration table to correct 
the difference in path lengths of two receivers in a dual¬ 
beam antenna with the ” horn-on-dish method” (Honma 
et al. 2008) (TBIN), and finally (10) setting the velocity 
axis (CVEL). Further details of the primary reduction are 
explained in Appendix 1 of Kurayama et al. (2011). 

After the primary reduction, we found three strong 
emission peaks detected in the cross-power spectrum, as 
shown in figure 2. These three peaks were found at the 
LSR velocity ranges of Vlsr = —32.31 to —32.89 km s _1 , 
—42.00 km s _1 , and —50.43 to —46.63 km s _1 in every 
observation. 

Since the typical synthesized beam-size of VERA is 
1 milli-arcsecond (mas) which is small compared to the 
known position accuracy of 15" (Sunada et al. 2007), we 
needed to make a fringe rate map (Walker 1981), which 
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enabled us to find the positions of maser spots with coarse 
resolution of ~ 0.1" in order to reveal the general distribu¬ 
tion of maser spots within the field of view. Note that a 
‘maser spot’ refers to a point like emission region detected 
in a particular velocity channel. Figure 3 shows the fringe 
rate map made using the channel of Vlsr = —32.31 km s -1 
during the first epoch. It shows that the position of the 
maser spots are shifted by Aacosi5 = —0.4 arcsecond from 
the tracking center in Right Ascension and AS = — 1.2 arc- 
second in Declination. The map nodes of other channels 
were almost the same as that of Vlsr = —32.31 km s -1 , 
within the coarse resolution. The fringe rate map did not 
change significantly throughout observations. 

Next we made dirty map data cubes using the AIPS 
task IMAGR to search for maser spots within a field 
of l."25 x 1"25 around the position determined with the 
fringe rate map, (Aa cos S, AS) = (—0"4, —1"2). We 
searched a velocity range of Vlsr = —56.33 to —23.04 
km s , to cover all the emitting components found in 
the range Vlsr = —50.43 to —32.31 km s -1 . We cata¬ 
loged maser spots found in each observation track and 
performed the CLEAN procedure (Hogborn 1974) for each 
maser spot using IMAGR. Finally, the position of each 
spot was measured using the AIPS task JMFIT. 

4. Results 

4-1. Detected, Maser Features 

Maser features were identified by the following criteria: 
the emission must be detected at a level of greater than 
five times the rms noise in individual channel maps and 
must appear over adjacent spectral channels. All detected 
maser features are listed in table 2, showing: (1) maser 
feature ID, (2) initial R.A. offset (AaocosiJ), (3) R.A. 
proper motion (/z a cos<5), (4) initial Dccl offset (A<5 0 ), (5) 
R.A. proper motion (ns), (6) LSR velocity (Vlsr), and 
(7) Detection epochs. 

The center of the maser distribution map was rede¬ 
fined as (a c ,S c ) = (21 h 39 m 40f501608, +51°20'32"60108) 
(J2000.0), which was the averaged position of the all de¬ 
tected maser features in the first epoch. 

Typically, more than 10 features were detected in most 
epochs, with the exceptions of epochs 5 and 8. In total, 
39 maser features were detected. Five features were iden¬ 
tified in six or more epochs and another 13 features were 
detected in two two four epochs. The other 21 features 
were detected only once. The identification criteria were 
(1) an identical maser feature moves less than 2.5 mas 
from its position in the previous observation, (2) a dif¬ 
ference in velocity giving peak temperature between con¬ 
secutive observations is less than 1 km s -1 , and (3) R.A. 
or Decl offsets must meet the theoretical curve of the an¬ 
nual parallax and linear proper motion within an error of 
1 mas. The first criterion excludes large proper motions 
exceeding 10 mas yr -1 considering a typical observation 
interval of 90 days. 


4-2. Annual Parallax and Proper Motion 

We used the five features (ID 1 -5 of table 2) detected 
in six or more epochs to calculate the annual parallax. 
Annual parallax fitting was performed for each maser 
spot. The annual parallax was calculated using 12 maser 
spots associated with 5 individual maser features. In ad¬ 
dition, we performed a combined-fitting in which the as¬ 
trometry of all maser spots are fit together assuming a 
common distance (see for example Sakai et al. (2012)). 
The combined-fitting produced a result which is of higher 
accuracy and consistent with that obtained for individual 
spots. The results are shown in table 3. Figure 4 shows 
the R.A. and Decl offsets of maser features plotted against 
time. The time variation of the R.A. offset Act cos (5^ and 
Decl offset A 5^ due to annual parallax and linear proper 
motion of ith feature is expressed as a function of time t , 
as follows. 

Aa cos S c ^ = w (— sin a c cos 0 + cose cos a c sin 0) ( 1) 

+/ 4 l) cos S c (t - 1^) + Aag ) cos^c 
AS (i) = w(— cosa c cos0 + cosesina c sin0)sin(5 c 
+ sine cos 5 C sin© 

+rf\t-tP) + A6P ( 2 ) 

where vo is the anuall parallax, e is the obliquity [e = 
23.440972), 0 is the ecliptic longitude of the Sun, 
and j/j 1 are proper motion in the R.A. and Decl direc¬ 
tions, is the first detection time of the ith feature, and 
Aaocosc)^ and AS^ are R.A. and Decl offsets at the time 
to- We fitted this theoretical curve to the observed posi¬ 
tions of the detected features by the least squares method, 
which is explained in Kovalevsky & Seidelmann (2012), 
with known parameters Aacos(5 c ^, A<?W, a c , S Cl 0, e, 
t, and t q \ to obtain unknown parameters w 1 p a cos<5 c ^\ 
(j!g\ AaoCOS(5 c ^ and ASq\ As seen in figure 4, most 
of the data were well fitted by the theoretical curve. We 
adopted positional errors of 0.24 mas, which were required 
to produce a reduced chi squared of one, and are indicated 
as the error bars in figure 4. 

The obtained annual parallax was 0.262 ± 0.031 mas, 
which corresponds to Heliocentric distance of D = 
3.82+Q4* kpc. Hereafter, we adopt this value as the dis¬ 
tance of IRAS 21379+5106. 

The averaged proper motions of 12 maser spots are 
also summarized in table 2. The proper motion of this 
system was calculated to be p a cosS = —2.74 + 0.08 mas 
yr _1 , ns = —2.87 + 0.18 mas yr _1 , which was the averaged 
value of the five maser features, adopting a standard er¬ 
ror. The systemic proper motion in Galactic coordinates 
was calculated to be /qcos6 = —3.95 + 0.18 mas yr _1 and 
pt = —0.34 + 0.08 mas yr _1 . Note that the error vector 
(A/j,icosb,Afj,b) = (0.18,0.08) mas yr _1 in the Galactic co¬ 
ordinate was calculated with the same coordinate trans¬ 
formation as the error vector {Afj, a cosS,Ans) = (0.08,0.18) 
mas yr” 1 in the equatorial coordinate. The corresponding 
linear velocity was Vi = Dpicosb = —71.5+ 10.1 km s _1 
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and Vb = Dfib = —6.21 ± 1.70 km s _1 in Galactic coordi¬ 
nates, when using D = 3.82 kpc. 

We also estimated proper motions of the other 13 maser 
features detected in only two to four epochs, in addition 
to the other five features, adopting the newly determined 
systemic distance of D = 3.82 kpc. 

Table 4 shows the internal motion of IDs 1-18, which 
were calculated by subtracting the systemic proper mo¬ 
tion, (fj, a cos5,]Zs) = (—2.74±0.08,—2.87±0.18) mas yr -1 
from the individual proper motions (/r Q cosd,/^), and by 
subtracting the systemic Vlsr = —42.27 ±0.22 km s -1 
from those of the individual masers. The result is dis¬ 
cussed in section 5.4. 

4-3. 3-D Motion Relative to the Sun in the LSR coordi¬ 

nate 

We calculated the three dimensional velocity vectors 
(U, V. W) with respect to the LSR, using the LSR velocity, 
obtained parallax, and proper motion, following Reid et 
al. (2009). As shown in figure 6, we defined the LSR veloc¬ 
ity coordinate so that U, V, and W axes were aligned with 
( l,b ) = (0,0), (90°,0), and (0,90°), respectively, though the 
U axis is opposite to the classical definition (Mihalas & 
Binney 1981). 

First, the LSR velocity Vlsr = —42.27±0.22 km s -1 was 
converted to the heliocentric velocity Vtiei = —56.43 ±0.22 
km s _1 with 

Vkei = Vlsr — (C^© td cos l ± V)| td sin l ) cos b — W|j td sin b (3) 

where (t/| td ,F| td , W| td ) is the standard solar motion de¬ 
fined as Uq A = 10.3 km s -1 , Vj td = 15.3 km s _1 , and 
Wf td = 7.7 km s-\ whose absolute value and direction 
are 20 km s -1 and a(1900) = 18 h , 6(1900) = ±30 d , respec¬ 
tively. 

Second, the linear velocity vector (V, V,) was converted 
to a 3-D vector, relative to the Sun, as follows. 

U — Uq — (Vriei cos b — I4sin6)cos/ — Vsinl (4) 

V — Vq = (Vneicosfe — Vb sin b) sin/ ± Vicosl (5) 

W — Wq — Vb cos b ± Vfiei sin b (6) 

The obtained 3-D velocity vector relative to the Sun in 
the LSR coordinate was U — Uq = 76.4 ±10.4 km s -1 , 
V — Vq = —49.7±1.0 km s _1 , and W-W© = -5.29 ±1.70 
km s -1 . 

5. Discussion 

5.1. Absolute 3-D Motion in the Galaxy and the Rotation 
Velocity 

We first discuss the 3-D motion of the maser source, 
IRAS 21379±5106, in the Galaxy, based on the obtained 
values of the parallax and proper motion. 

5.1.1. 3-D Vector in the LSR coordinate 

The 3-D motion of the object in the LSR coordi¬ 
nate was calculated by adding the solar peculiar motion 
(Uq,Vq,Wq) to the relative velocity obtained in the pre¬ 
vious section. 


Thanks to the recent improvements that come with 
the continuous developments in observational techniques 
and methods, the solar peculiar motion relative to the 
LSR has been measured more accurately than the classi¬ 
cal value of (C/,| td , Vj td , W<| td ). Dehnen & Binney (1998) 
estimated it to be Uq = 10.0 km s -1 , Vq = 5.2 km s _1 , 
and Wq = 7.2 km s -1 based on the parallax and proper 
motions obtained from the Hipparcos catalog (Perryman 
& ESA 1997). Schonrich et al. (2010) estimated the so¬ 
lar 3-D motion of the LSR so that it should be insen¬ 
sitive to the metallicity gradient and concluded it to be 
Uq = 11.1 km s -1 , Vq = 12.2 km s -1 , and Wq = 7.3 km 
s _1 . Honma et al. (2012) took the same values for Uq and 
Vq as Dehnen & Binney (1998) and measured Vq to be 
12.0 km s _1 by analyzing VLBI astrometry data for 52 
Galactic maser sources. Comparing these 3 papers with 
each other, it can be regarded that Uq and Wq are well 
determined. Regarding the other component Vq, the re¬ 
sults from Schonrich et al. (2010) and Honma et al. (2012) 
are consistent with each other. Therefore, we adopt the 
solar peculiar motion relative to the LSR as Uq = 10.0 km 
s _1 , Vq = 12.0 km s -1 , and Wq = 7.2 km s _1 as adopted 
by Honma et al. (2012). We conclude that the 3-D veloc¬ 
ity vector of the IRAS 21379±5106 in the LSR coordinate 
is U = 86.4 ± 10.1 km s -1 , V = —37.7 ± 1.0 km s -1 , and 
W = 1.91 ± 1.70 km s" 1 . 

5.1.2. Absolute 3-D Motion in the Galaxy 

With the 3-D vector of the source in the LSR coor¬ 
dinate, we can estimate its absolute 3-D motion in the 
Galaxy. 

First, the Galactocentric distance to the source is cal¬ 
culated using the following equation, 

R = Rq + ( Dcosb ) 2 — 2RoDcosbcosl. (7) 

Next we define the velocity vector as (Vft in , Vg, V z ), 
where V/j in ,Vg,I4 are the velocity components directed 
to the Galactic center, the rotational velocity about the 
Galactic center, and the vertical velocity perpendicular to 
the Galactic plane (the direction of the north Galactic pole 
is positive), respectively, as outlined in figure 6. Using this 
definition, each component is estimated using the follow¬ 
ing equations, 

VR ln =f7cos0- (V±Vo)sin<?i (8) 

Vg = (V ± Vo)cos<)>± U sin<)> (9) 

V, = w (10) 

where </> is the angle between the Sun and the source as 
viewed from the Galactic center given by 

. , Dcosb R 0 — Dcosbcosl , 

sm0=—-—sm l cos cp= -—-. (11) 

R R 

Since the velocity vector (VR in ,Vg,T4) depends on 
the Galactic constants, Rq and Vo, which are the 
Galactocentric distance of the Sun and the rotational ve¬ 
locity of the Sun about the Galactic center, respectively, 
we estimated the 3-D velocity vector using some possible 
values reported in recent studies. 

First we took the Galactic constants from Honma et 
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al. (2012), who concluded that the Galactic constants are 
Ro = 8.05 kpc and Vo = 238 km s _1 by analyzing the avail¬ 
able latest astrometric data obtained with VLBI includ¬ 
ing VERA, VLBA, and EVN. Using these Galactic con¬ 
stants we estimate the Galactocentric distance of IRAS 
21379+5106 to be R = 9.22 kpc and the 3-D velocity 
components, Vn in = —3.8 km s -1 , Vq = 218 km s -1 , and 
14 = 1.91 km s -1 . 

If we take IAU recommendation value, Rq = 8.5 kpc 
and Vo = 220 km s” 1 (Kerr & Lynden-Bell 1986), which 
are the most commonly adopted values, we estimate a 
Galactocentric distance of R = 9.63 kpc and the velocity 
vector components Vn in = 7.5 km s -1 , Vq = 202 km s” 1 , 
and V z = 1.91 km s -1 . 

According to Dehnen & Binney (1998), a mass model 
with Ro = 8.0 kpc and Vo = 217 km s -1 is best fitted to the 
observational data. Adopting these Galactic constants, we 
estimate a Galactocentric distance of R = 9.18 kpc and 
velocity vector components V/j in = 4.4 km s _1 , Vq = 199 
km s _1 , and V z = 1.91 km s . 

Finally when we adopt the Galactic constants Ro = 8.34 
kpc and Vo = 240 km s -1 , recently estimated by Reid 
et al. (2014) by analyzing astrometry data with VLBA, 
the Galactocentric distance and the velocity vectors are 
estimated as R = 9.48 kpc, V/j in = —1.9 km s -1 , Vq = 220 
km s -1 , and V z = 1.91 km s -1 , respectively. These results 
are summarized in tabic 5. 

These results systematically show that the rotational 
velocity Vq tends to be slower than that of the Sun. This 
is consistent with the declining rotation curve beyond 
R = Ro obtained with the recent work on maser sources 
associated with the Perseus arm (Sakai et al. 2012; Oh 
et al. 2010) and is consistent with the existence of 9- 
kpc dip found in the Galactic rotation curve (Sofue et al. 
2009), which might be caused by the mass excess at the 
Galactocentric distance of 11 kpc. Regarding the radial 
velocity component, VA in , the error was not small enough 
to judge which direction the object moves. The vertical 
component V z = 1.91 km s _1 is small enough to be con¬ 
sistent with the natural random velocity, whose typical 
velocity dispersion is 2 11 km s -1 (Malhotra 1994). 

5.2. Spectral Type of YSO 

Our accurately measured distance modifies the bolo- 
metric luminosity of the central young stellar object 
(YSO). Wang et al. (2009) estimated it to be 5.18 x 10 3 Lq, 
which corresponds to the spectral type B1 based on 
IRAS data and a kinematic distance of 5.9 kpc, when 
it is assumed to be a Zero Age Main Sequence star 
(ZAMS)(Panagia 1973). With the newly determined 
distance, the bolometric luminosity was modified to be 
(2.15 + 0.54) x 10 3 Lq by multiplying with the square of 
the distance ratio (3.8/5.9) 2 . As a result, we re-estimate 
the spectral type of the central YSO to be B2-B3. This 
result is still consistent with the fact that this object is 
known to have an ultra compact H II region, such objects 
are associated with YSOs of main sequence luminosity 
corresponding to type B3, or hotter. This picture is valid 
even though the object was found to be less luminous than 


previously thought. 

5.3. Distance to the Perseus Arm 


Since IRAS 21379+5106 is a member of the Perseus arm 
of the Milky Way Galaxy, we also consider that its trigono¬ 
metric distance can be regarded as that of the Perseus 
arm. We compare our result with other relevant distance 
values. 

1. Another VLBI astrometric measurement: The annual 
parallax of IRAS 21379+5106 was recently reported by 
Choi et al. (2014), who reported an annual parallax of 
0.206 + 0.007 mas and a corresponding trigonometric dis¬ 
tance of 4.85to'i6 kpc, b ase d on astrometric observations 
with the VLBA. The difference in parallax values is about 
2cr. A possible reason for this discrepancy is that they 
used only a single feature rather than the cumulative re¬ 
sults of multiple maser features, as was done in this paper. 

2. Kinematic Distance: The kinematic distance is often 
used as a distance estimation and Wang et al. (2009) cal¬ 
culated it to be 5.9 kpc. If we assume a flat rotation 
curve, the Galactocentric distance R and the Heliocentric 
distance D are calculated using the LSR velocity as fol¬ 
lows. 


R — Rq 


_Vb_1_ 

Vlsr sin/cos 6 



D = RqcosI + 



Rq sin 2 1 


( 12 ) 

(13) 


When we adopted Galactic constants we use in this paper, 
Ro = 8.05 ± 0.45 kpc and Vo = 238 ± 14 km s _1 , we have 
R = 9.78 ± 0.56 kpc and D = 4.86 ± 1.30 kpc, which turn 
out be overestimations, compared to our trigonometric 
distance. Adopting any of the other Galactic constants 
mentioned in the previous section, the kinematic distance 
is still larger than the trigonometric distance, as long as 
we assume a flat rotation curve. 

3. S124'- Sharpless 124 (S124) is an H n region seen in 
figure 1 as a circular-shaped emission region just below 
IRAS 21379+5106. With celestial coordinates and LSR 
velocity (/,&) = (94.°57, —1.°45) and Vlsr = —43.4 km s _1 , 
it is closely associated with IRAS 21379+5106, thus these 
regions are thought to coexist. Brand & Blitz (1993) 
estimated its distance to be 2.6 ± 0.6 kpc by photomet¬ 
ric distance, which underestimated it by 30%. However, 
our newly determined trigonometric distance seems more 
plausible because applying the photometric distance of 
S124 to IRAS 21379+5106 would re-classify the central 
YSO as a spectral type later than B3, which would not 
be hot enough to produce the observed ultra compact H 
II region. 

4■ AFGL 2789: AFGL 2789 is another object located 
in the Perseus arm, and its position (/,&) = (94.°6, —1.°80) 
and LSR velocity (—44 km s _1 ) are close to those of IRAS 
21379+5106 (Oh et al. 2010). This object is also seen in 
figure 1 as a bright point source located at the south¬ 
ern edge of the circular-shaped infrared emission of S124. 
Oh et al. (2010) concluded that trigonometric distance of 
AFGL 2789 was 3.08 + 0.30 kpc, which is consistent with 
our study, within the error. 
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5-4- Asymmetric Spectrum and Maser Distribution - 
Evidence of Bipolar Outflow 

The maser distribution shown in figure 5 indicates that 
most of the maser features seem to be linearly aligned 
on the sky plane. North-east and south-west components 
tend to be red- and blue-shifted, respectively. We drew a 
straight line in figure 5, which is fitted with the positions 
of all the maser features listed in table 4 except IDs 3, 4, 
8, and 9 which deviate largely from the line. The observed 
situation resembles a bipolar outflow from the YSO, which 
is reasonable considering this object is an ultra-compact 
H II region and high-mass star-forming region (Bronfman 
et al. 1996; Wang et al. 2009; Wu et al. 2006; Shepherd & 
Churchwell 1996), which are often associated with molec¬ 
ular outflows (Ortega et al. 2012; Qin et al. 2008; Tofani 
et al. 1995). 

The blue-shifted component is dominant in the maser 
distribution map. This corresponds to the asymmetry of 
the spectrum shown in figure 2, where the blue compo¬ 
nent is brighter than the red one. This asymmetry is also 
found in the 12 CO(J =1-0) and 13 CO(J =3-2) spectra 
obtained with the IRAM and KOSMA telescopes, respec¬ 
tively (Wouterloot & Brand 1989; Wang et al. 2009). 

A plausible explanation of this characteristics is that 
there exists a disk around the central young star and the 
red component is obscured by this optically thick disk, as 
seen in Motogi et al. (2013). 

Next let us consider the internal motions of the maser 
features shown in table 4 and figure 7. Maser features 
seem to be ejected from the central YSO, which is thought 
to be located on the fitted line. These are the similar kine¬ 
matic characteristics as were found in the proper motion 
of water-vapor masers associated with IRAS 21391+5801 
(Patel et al. 2000). In a further detail, the directions of 
vectors are not always aligned with the fitted line. This 
implies that the motions of the maser features are not 
simple ballistic motions from the central YSO, but seem 
perturbed possibly due to interaction with ambient mat¬ 
ter (Lada & Fich 1996; Yu et al. 1999) and/or helical 
motion driven by a magnetic field (Matzner & McKee 
1999; Woitas et al. 2005). 

The maser features located in the region of Decl offset 
— 100 mas to —150 mas seem to move eastwards, differ¬ 
ently from the bipolar flow motion traced by the other 
maser features. A possible explanation for this discon¬ 
tinuity is that the bipolar flow hits against the ambient 
interstellar medium which subsequently changed the di¬ 
rection of the outflow. 

6. Summary 

We measured the annual parallax and proper motion of 
the high-mass star-forming region IRAS 21379+5106 with 
VERA. The resultant annual parallax was 0.262 ± 0.031 
mas, which corresponds to a trigonometric distance of 
3.82j/o;4i kpc. The proper motion was p a cos5 = —2.74 + 
0.08 mas yr -1 and ps = —2.87 + 0.18 mas yr -1 . When 
we adopt the Galactic constants i?o =8.05 + 0.45 kpc and 
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Vo = 238+ 14 km s -1 , suggested by Honma et al. (2012), 
we conclude that the Galactocentric distance of IRAS 
21379+5106 is R = 9.22 ± 0.43 kpc and the velocity com¬ 
ponents in cylindrical coordinates are V/j in = —3.8 + 16.9 
km s -1 , Vg = 218 + 19 km s -1 , and V z = 1.91 + 1.70 km s -1 . 
Using our new distance, the spectral type of the central 
YSO was re-estimated to be B2-B3, but not Bl, which was 
its previous classification. Masers are distributed along a 
straight line indicating the existence of a bipolar outflow. 
All the derived physical parameters are summarized in 
table 6. 

We thank the VERA project members for supporting 
observations. We would also like to thank the referee for 
carefully reading the manuscript. The authors acknowl¬ 
edge financial support from JSPS/MEXT KAKENHI 
Grant Numbers 26800104, 13J03569, 25610043, 24540242, 
and 25120007. 
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Fig. 1. Mid-Far infrared composite image of 7.° 5 x 7.° 5 area 
centered on IRAS 21379+5106. This was made using IRIS 25 
//m-, 60 //m-, and 100 |zm-band images. IRAS 21379+5106 is 
indicated with a green circle. A circular-shaped emission re¬ 
gion just below IRAS 21379+5106 is the H II region Sharpless 
124 (S124). A bright point source located at the southern 
edge of the circular-shaped infrared emission of SI24 is AFGL 
2789. 
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Fig. 2. Cross-power spectrum obtained at the baseline be¬ 
tween Mizusawa and Iriki at UT 7:39:47.3 of 6th November 
2009. The channel range is from 140 (Vlsr = —14.19 km 
s _1 )to 260 (Vlsr = —64.75 km s _1 ). 
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Fig. 4. Top: Time variations of R.A. (left) and Decl (right) offsets. The potted curves are the best-fitted models including annual 
parallax and linear proper motion. Bottom: Time variation of R.A. (left) and Decl (right) offsets with the proper motions subtracted. 
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Fig. 5. Map of maser features listed in table 2. The color denotes the line-of-sight velocity and its value is indicated in the vertical 
color bar. Each symbol denotes the observation epoch shown at the top-right corner. Each feature shows movement between epochs 
due to the proper motion and annual parallax. The angular size of 1 mas corresponds to a physical size of 3.82 AU. 



Fig. 6. Schematic depiction definining the velocity vectors (U,V) and (Vr^ t Vg). 


Table 1. Summary of Observations 


Epoch 

Date 

DOYT 

Synthesized beam 
(mas x mas) 

PA 

(deg) 

Noise 
(mJy beam -1 ) 

1 

2009 NOV 06 

310 

1.21 x 0.75 

-58 

7.3 

2 

2010 JAN 21 

386 

1.18 x 0.75 

-48 

8.7 

3 

2010 MAR 15 

439 

1.20 x 0.77 

-60 

13.3 

4 

2010 MAY 28 

513 

1.16 x 0.77 

-48 

9.7 

5 

2010 OCT 13 

651 

1.21 x 0.73 

-41 

3.8 

6 

2010 DEC 09 

708 

1.21 x 0.74 

-41 

4.2 

7 

2011 FEB 27 

788 

1.13 x 0.73 

-48 

6.3 

8 

2011 AUG 02 

944 

1.14 x 0.74 

-43 

13.2 

9 

2011 NOV 05 

1039 

1.14 x 0.73 

-39 

14.5 


t DOY (Day-of-year) indicates days from 2009 January 1. 
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Table 2. List of Maser features 


ID 

Aao cos S 

/To, cos S 

A 6 o 

H5 

Vlsr 

Detection 


mas 

mas yr -1 

mas 

mas yr -1 

km s " 1 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 

(7) 

1 

74.228 ±0.213 

-2.226 ±0.141 

206.406 ±0.212 

-2.201 ±0.141 

-32.3 

123456*** 

2 

71.831 ±0.155 

-2.866 ±0.083 

194.856 ±0.156 

-3.317 ±0.083 

-35.7 

123456789 

3 

22.334 ±0.155 

-2.987 ±0.083 

-131.761 ±0.156 

-3.211 ±0.083 

-42.0 

123456789 

4 

40.727 ±0.155 

-3.260 ±0.083 

-127.782 ±0.156 

-3.703 ±0.083 

-42.8 

123456789 

5 

-57.609 ±0.470 

-4.389 ±0.411 

8.706 ±0.457 

-5.244 ±0.398 

-43.7 

*234567** 

6 

-10.550 ±0.470 

-2.726 ±0.411 

47.374 ±0.457 

-1.793 ±0.398 

-47.9 

1234***** 

7 

-10.283 ±0.273 

-2.849 ±0.166 

44.553 ±0.272 

-2.189 ±0.166 

-49.2 

1234***** 

8 

-20.904 ±0.070 

-2.683 ±0.973 

-144.552 ±0.058 

-3.102 ±0.495 

-45.4 

12 ******* 

9 

-18.711 ±0.126 

-0.122 ±0.891 

-144.697 ±0.050 

-3.045 ±0.372 

-46.2 

^ 2 ******* 

10 

71.756 ±0.053 

-4.402 ±1.207 

200.915 ±0.053 

-3.179 ±0.679 

-32.3 

*23****** 

11 

-13.830 ±0.136 

-1.273 ±0.691 

37.586 ±0.148 

-2.598 ±0.743 

-46.6 

**34***** 

12 

-97.535 ±0.051 

-1.040 ±0.373 

-128.411 ±0.065 

-6.251 ±0.454 

-39.1 

****^ 0 *** 

13 

-48.094 ±0.050 

-2.858 ±0.298 

-17.912 ±0.053 

-8.050 ±0.293 

-32.7 

***** 07 ** 

14 

-96.495 ± 0.034 

-6.743 ±2.735 

-132.491 ±0.037 

6.003 ±3.173 

-40.3 

***** 07 ** 

15 

-16.759 ±0.539 

-3.202 ±1.276 

26.884 ±0.534 

-3.398 ±1.262 

-44.9 

******7g* 

16 

-39.055 ±0.152 

-3.114 ±0.709 

24.038 ±0.064 

-7.305 ±0.745 

-50.0 

***** 07 ** 

17 

-108.001 ±0.207 

-0.116 ±1.042 

-118.739 ±0.270 

-5.270 ±1.407 

-39.1 

*******gg 

18 

-45.574 ±0.185 

-1.696 ±0.783 

24.077 ±0.169 

-4.851 ±0.711 

-49.2 

*******gg 

19 

-55.938 ±0.068 


27.719 ±0.055 


-41.2 

1 ******** 

20 

-13.798 ±0.076 


40.478 ±0.062 


-48.7 

i******** 

21 

-13.833 ±0.130 


40.094 ±0.071 


-49.6 

j******** 

22 

-21.679 ±0.155 


42.460 ±0.070 


-44.9 

* 2 ******* 

23 

-39.733 ±0.124 


-28.275 ±0.061 


-43.7 

***^***** 

24 

-36.587 ±0.033 


27.239 ±0.022 


-50.0 

***^ 2 ***** 

25 

-57.389 ±0.030 


19.943 ±0.032 


-34.0 

***** 0 *** 

26 

69.171 ±0.085 


193.120 ±0.106 


-36.9 

***** 0 *** 

27 

20.905 ±0.036 


-152.029 ±0.038 


-41.2 

***** 0 *** 

28 

-15.774 ±0.057 


33.159 ±0.126 


-46.6 

***** 0 *** 

29 

-18.188 ±0.041 


26.321 ±0.089 


-47.5 

********g 

30 

-41.598 ±0.068 


77.481 ±0.052 


-25.6 

******7** 

31 

69.312 ±0.100 


204.313 ±0.046 


-29.4 

******7** 

32 

-16.799 ±0.146 


26.571 ±0.064 


-47.1 

******7** 

33 

-16.135 ±0.096 


39.354 ±0.187 


-48.7 

******7** 

34 

-16.316 ±0.302 


27.892 ±0.395 


-53.0 

******7** 

35 

-56.108 ±0.170 


23.175 ±0.138 


-37.8 

>|ofc*****g* 

36 

-44.129 ±0.067 


70.453 ±0.088 


-46.6 

>|ofc*****g* 

37 

-42.198 ±0.107 


18.659 ±0.202 


-50.0 

>|ofc*****g* 

38 

-99.698 ±0.069 


-136.305 ±0.033 


-41.2 

********g 

39 

—18.684±0.110 


32.914 ±0.089 


-47.1 

********g 


(1) maser feature ID, (2) initial R.A. offset (Aaocos<5) from a c = 21 h 39 m 40. s 501608, (3) R.A. proper 
motion (/x a cos 8), (4) initial Decl offset (ASo) from 5 C = ±51°20'32"60108, (5) R.A. proper motion (ps), 
( 6 ) LSR velocity (Vlsr), and (7) Detection epochs, where the symbol * indicates epochs the maser 
feature was not detected. Since maser features ID 19-39 were detected only once, the proper motion in 
colums (3) and (4) are not calculated. 
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Fig. 7. Vector map of the internal motion made with the proper motions listed in table 4. The internal motion is obtained by 
subtracting the averaged motion individual feature motions. The black thick arrow at the bottom-left indicates vector of 100 km 
s _1 . The angular size of 1 mas corresponds to a physical size of 3.82 AU. 


Table 3. Parallax and Proper Motion 


ID 

Vlsr 

Parallax 

Aao cos S 


Ha COS S 



AS 0 



Hs 


(km s _1 ) 


i 

(mas) 



(mas) 

(mas 

yr !) 



(mas) 

(mas 

yr- 1 ) 

(1) 

(2) 



(3) 



(4) 



(5) 



(6) 



(7) 

la 

-31.5 

0.201 

± 

0.152 

74.392 

± 

0.478 

-2.377 

± 

0.356 

206.438 

± 

0.480 

-2.048 

± 

0.356 

lb 

-31.9 

0.256 

± 

0.129 

74.310 

± 

0.407 

-2.308 

± 

0.303 

206.295 

± 

0.408 

-1.993 

± 

0.303 

lc 

-32.3 

0.318 

± 

0.108 

74.212 

± 

0.286 

-2.237 

± 

0.189 

206.069 

± 

0.283 

-2.016 

± 

0.190 

2a 

-34.8 

0.222 

± 

0.172 

71.284 

± 

0.802 

-2.668 

± 

0.345 

195.872 

± 

0.811 

-3.656 

± 

0.354 

2b 

-35.3 

0.243 

± 

0.092 

71.839 

± 

0.229 

-2.880 

± 

0.123 

194.772 

± 

0.239 

-3.241 

± 

0.126 

2c 

-35.7 

0.264 

± 

0.091 

71.823 

± 

0.226 

-2.842 

± 

0.121 

194.807 

± 

0.235 

-3.355 

± 

0.124 

2 d 

-36.1 

0.257 

± 

0.111 

71.876 

± 

0.304 

-2.844 

± 

0.157 

194.799 

± 

0.320 

-3.381 

± 

0.161 

3a 

-41.6 

0.381 

± 

0.091 

22.263 

± 

0.270 

-2.903 

± 

0.174 

-132.067 

± 

0.270 

-2.953 

± 

0.174 

3b 

-42.0 

0.263 

± 

0.097 

22.315 

± 

0.240 

-2.980 

± 

0.128 

-131.789 

± 

0.250 

-3.201 

± 

0.132 

4a 

-42.8 

0.274 

± 

0.100 

40.686 

± 

0.249 

-3.242 

± 

0.133 

-127.820 

± 

0.259 

-3.692 

± 

0.137 

5a 

-48.7 

0.196 

± 

0.119 

-10.302 

± 

0.414 

-2.838 

± 

0.253 

44.636 

± 

0.409 

-2.325 

± 

0.246 

5b 

-49.2 

0.263 

± 

0.111 

-10.358 

± 

0.387 

-2.793 

± 

0.237 

44.948 

± 

0.383 

-2.495 

± 

0.231 


Combined 0.262 ±0.031 -2.741 ±0.082 -2.865 ±0.180 


(1) maser spot ID, (2) LSR velocity, (3) obtained parallax, (4) initial R.A. offset (Aaocos<5) from a c = 21 h 39 m 40?501608, 
(5) R.A. proper motion (/z a cos 5), (6) initial Decl offset (A<5o) from 5 C = +51°20'32"60108, and (7) R.A. proper motion 
(hs), (6) radial velocity (V)- The last row shows the result of combined-fitting, which is the annual parallax fitting 
made for all maser spots simultaneously assuming the common distance. 
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Table 4. Internal Motions 


ID 

R.A. 

Ha COS S c — fi a COS S c 

(mas yr -1 ) 

v a -v a 

(km s -1 ) 

Decl 

M(5 - JjJ Vs - Vs 

(mas yr -1 ) (km s -1 ) 

Vlsr 

(km s -1 ) 

1 

0.52 ±0.27 

9.3 ±5.0 

0.66 ±0.23 

12.0 ± 4.1 

-32.3 ±0.4 

2 

-0.13 ±0.30 

—2.3 ± 5.4 

-0.45 ±0.20 

-8.1 ±3.7 

—35.7 ± 0.4 

3 

-0.25 ±0.30 

-4.5 ±5.4 

-0.34 ±0.21 

-6.2 ±3.7 

-42.0 ±0.4 

4 

-0.52 ±0.32 

—9.4 ±5.7 

-0.84 ±0.20 

-15.1 ±3.6 

-42.8 ±0.4 

5 

-1.64 ±0.48 

-29.5 ± 8.6 

-2.38 ±0.21 

-42.9 ±3.8 

-43.7 ±0.4 

6 

0.03 ±0.56 

0.5 ±10.2 

1.07 ± 1.54 

19.4 ±27.9 

-47.9 ±0.4 

7 

-0.11 ±0.41 

—2.0 ± 7.3 

0.68 ±0.30 

12.2 ±5.4 

-49.2 ±0.4 

8 

0.07 ±0.98 

1.2 ± 17.6 

-0.25 ±0.53 

-4.6 ±9.5 

-45.4 ±0.4 

9 

2.63 ±0.89 

47.5 ±16.2 

-0.20 ±0.41 

-3.5 ±7.5 

-46.2 ±0.4 

10 

-1.64 ±1.21 

—29.6 ± 21.9 

-0.31 ±0.70 

—5.6 ± 12.7 

-32.3 ±0.4 

11 

1.47 ±0.70 

26.6 ± 12.6 

0.29 ±0.76 

5.2 ±13.8 

-46.6 ±0.4 

12 

1.70 ±0.38 

30.7 ± 6.9 

-3.41 ±0.49 

-61.5 ± 8.8 

-39.1 ±0.4 

13 

-0.10 ±0.31 

— 1.8 ± 5.6 

-5.19 ±0.34 

-93.7 ±6.2 

-32.7 ±0.4 

14 

-3.98 ±2.74 

-72.0 ±49.4 

8.86 ±3.18 

160.1 ±57.4 

-40.3 ±0.4 

15 

-0.46 ±1.28 

-8.3 ±23.1 

-0.52 ±1.27 

-9.4 ±23.0 

-44.9 ±0.4 

16 

-0.35 ±0.71 

—6.4± 12.9 

-4.45 ±0.77 

-80.3 ±13.8 

-50.0 ±0.4 

17 

2.61 ±1.05 

47.1 ±18.9 

-2.42 ±1.42 

-43.7 ±25.6 

-39.1 ±0.4 

18 

1.03 ±0.79 

18.6 ±14.2 

-2.00 ±0.73 

-36.1 ±13.2 

-49.2 ±0.4 


Table 5. 3-D velocity in the Cylindrical Coordinate Adopting Some Possible Galactic Constants 


Ro 

kpc 

Vo 

km s -1 

R 

kpc 

v Rin 

km s 1 

V e 

km s 1 

v z 

km s 1 

Reference 

8.05 ±0.45 

238 ± 14 

9.22 ±0.43 

-3.8 ±16.9 

218 ±19 

1.91 ±1.70 

Honma et al. (2012) 

8.5 

220 

9.63 ±0.13 

7.5 ±13.5 

202 ±7 

1.91 ±1.70 

Kerr & Lynden-Bell (1986) 

8.0 

217 

9.18 ± 0.17 

4.4 ± 13.7 

199 ±7 

1.91 ±1.70 

Dehnen & Binney (1998) 

8.34 ±0.16 

240 ±8 

9.49 ±0.20 

-1.9 ±14.9 

220 ±11 

1.91 ±1.70 

Reid et al. (2014) 


Table 6. Summary of Adopted and Derived Physical Parameters 


Phase Tracking Center 

(a, S) = (21 h 39 m 40f55, ±51°20'34"0) (J2000.0) 

(l,b) = (95.°29670, —0.°93660) 

Center of Maser Map 

(a c ,S c ) = (21 h 39 m 40f501608, ±51°20'32"60108) (J2000.0) 
(/ C ,6 C ) = (95.°29635672, —0.°93695671) 

LSR Velocity 

Vlsr = -42.27 ± 0.22 km s" 1 

Heliocentric Velocity 

V Hel = -56.43 ± 0.22 km s” 1 

Annual Parallax 

vo = 0.262 ± 0.031 mas 

Trigonometric Distance from the Sun 

D = 3.82±g-|f kpc 

The Galactic constants 

Ro = 8.05 ± 0.45 kpc, V 0 = 238 ± 14 km s"^ 

3-D Solar Motion in LSR 

(U 0 ,V Q ,W Q ) = (10.0,12.0,7.2) km s" 1 * 

Galactocentric Distance 

R = 9.22 ± 0.43 kpc 

Proper Motion 

(fx a cos S,fj,s) = (—2.74 ±0.08,—2.87 ±0.18) mas yr -1 
(/x/cosfc,^;,) = (—3.95 ± 0.18, —0.34 ± 0.08) mas yr -1 
(Vi,V b ) = (-71.5 ±10.1,-6.21 ±1.70) km s” 1 

3-D Velocity Vector 

( U,V,W) = (86.4± 10.1,-37.7± 1.0,1.92± 1.70) km s” 1 
lv Rin ,V e ,V z ) = (—3.8± 16.9,218± 19,1.92± 1.70) km s” 1 

Bolometric Luminosity of the YSO 

(2.15 ±0.54) x 10 3 L q 

Spectral Type of the YSO 

B2-B3 


T Honma et al. (2012) 
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Fig. 3. Fringe-rate map of the maser emission detected at 
the channel of Vlsr = —32.31 km s —1 in the first epoch. 
The map origin is located at the phase tracking center 
(a, 5 ) = (21 : 39 : 40.55,+51: 20 : 34.0) (J2000.0). The position 
of the maser spot is shifted by Aa:cos<5 = —0.4 mas from the 
phase tracking center in Right Ascension and A5 = —1.2 mas 
in Declination. 



